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So far, the infrared behavior of the gluon and ghost propagator based on the Gribov-Zwanziger approach pre- 
dicted a positivity violating gluon propagator vanishing at zero momentum, and an infrared enhanced ghost 
propagator. However, recent data based on huge lattices have revealed a positivity violating gluon propagator 
which turns out to attain a finite nonvanishing value very close to zero momentum. At the same time the ghost 
propagator does not seem to be infrared enhanced anymore. We point out that these new features can be ac- 
counted for by yet unexploited dynamical effects within the Gribov-Zwanziger approach, leading to an infrared 
behavior in qualitatively good agreement with the new data. 

PACS numbers: ll.lO.Gh, 12.38.Aw, 12.38.Lg 



I. INTRODUCTION. 

During the past few years, much attention has been devoted 
to the study of the gluon and ghost propagator, including 
their low energy behavior where Yang-Mills gauge theories 
are confining. As a consequence the gluon cannot be consid- 
ered as a free particle anymore. Due to the lack of an explicit 
knowledge of the physical degrees of freedom, it still remains 
highly useful to study the gluon and ghost propagator in or- 
der to probe the nonperturbative infrared regime. Let us only 
mention that the gluon propagator for example finds an impor- 
tant use in phenomenological studies, see e.g. U]]. In the past, 
a good agreement between the lattice data and the analytical 
results arising from the Gribov-Zwanziger action in the Lan- 
dau gauge were found: (1) an infrared suppressed and positiv- 
ity violating gluon propagator vanishing at zero momentum, 
(2) an infrared enhanced ghost propagator. We recall that the 
Gribov-Zwanziger action was constructed to take into account 
the existence of gauge copies 12^ JJ. However, very recent lat- 
tice data obtained at large volumes fl \^, which allows one 
to get very close to zero momentum, now give evidence of 
a hitherto unexpected behavior in the deep infrared: (1) an 
infrared suppressed and positivity violating gluon propagator 
nonvanishing at zero momentum, (2) a ghost propagator essen- 
tially behaving like 1 / p^ at low momentum, which is clearly 
not enhanced. To our knowledge, none of the current theo- 
retical approaches exhibit all such features la, Q H B llOtl . 
In this letter we propose a dynamical mechanism within the 
Gribov-Zwanziger approach that could account for the new 
lattice results. 
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II. THE GRIBOV-ZWANZIGER ACTION. 

We first give a short overview of the action constructed by 
Zwanziger |3] which implements the restriction to the Gri- 
bov region £1 in Euclidean Yang-Mills theories in the Landau 
gauge. We recall that this restriction to £1 can be implemented 
by adding the nonlocal horizon function to the original Yang- 
Mills action. 



5yM + ^Landau " f g" f d^xf^'^'A^ {M -')"' /"A^ , 



(1) 



where M"" = -3^, (3^5"* + g/"'"A;,) is the Faddeev- 

Popov operator, 5ym = 1 /4 / d'^xfuv-F^v and ^Landau = 
/ A'^x{b"dijA" + c"df^fl"j'c'') stands for the gauge fixing and the 
ghost part. However, it is unclear how to handle consistently 
such a nonlocal action at the quantum level, so we are obliged 

to add extra fields ((p^'',(p][;^',co^'',(0^' j in order to localize this 

action. Doing so, the Gribov-Zwanziger action reads IBI UOII 

s=so-y'gld'x (^r'"A>;r+r^^A;,(p^'+l (A^2_ j)^^ ^ 

(2) 
with 



So 



Sym + 1 d^x {b"d,Al + rd, (D,cf) 



(3) 

whereby (^'^,92'^ ) are a pair of complex conjugate bosonic 
fields, whereas (co"*^^, co"'^) are anticommuting. Due to a global 
U (/) symmetry, / = 4 (A^^ ~ 1 ) > with respect to the composite 



index / = (/i, c) of the additional fields ( cp^ 
we introduced a notational shorthand. 



^'p;;' 



CO r , CO 



,W,<9";,^, 
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((p«,(p?,cof,cof) 



(4) 



The dimensional parameter y is not free, being determined by 
the gap equation (horizon condition) dF/dy = 0, which en- 
sures the restriction to the Gribov region. F is the quantum 
effective action obtained from (|2]l. 

As it has been shown in [3^, J^, the action dU is renormaliz- 
able to all orders. To prove this by the method of algebraic 
renormalization [11], we embed this action into a larger one 
which contains more symmetries. This action turns out to be 
given by ifToll 

E = So + Ss + Sen, (5) 

with ^0 given in (O and 



5ext = / d'x (~K"^ [D^cf + ^-gL^f'^c'A . (6) 



We introduced new sources M^", V^' ,f/^", N"^' , K", and L", 
which are necessary to analyze the renormalization of the 
corresponding composite field operators in a BRST invariant 
fashion. The BRST operator s acts on the fields and sources 
appearing in the action as follows 



sA"^ = - iD,cY' ,sc" = -gf^c'c' ,sc" ^b\sb"^0, 



si^ = coj' ,.cof = ,smi = cpf ,^cpf - ,su;' 



m: 



sM" 



: , sVf = N'^' , sN'^' = , sK^ = 0, sL" = , 



(7) 



whereby the BRST operator s is nilpotent, s-^ — 0. One can 
easily see that the action E is indeed BRST invariant. Hence- 
forth, the action E displays a greater number of symmetries, 
encoded in the following Ward identities JSl [l^l ■ 



For the U (/) invariance mentioned before we have 



UijT. = 0, 



(8) 



Ui, = 



I'''(^'wr^'^^^'ii-^'i^] 



The Slavnov-Taylor identity reads 

5(E) =0, 



(9) 






• The Landau gauge condition and the antighost equation 
are given by 



^ -dA" 
db" ~ " '" 

5E ^ 5E 



(10) 
(11) 



• The ghost Ward identity: 



(12) 
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A^i = g fd*xf^'(K'^A'^~d'c' 



The term A"[ denotes a classical breaking as it is linear 
in the quantum fields. 

• The linearly broken local constraints: 
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(15) 
(16) 



8E 



-8f""'u';-^ = gf'XK+J^! 



• The exact ^J^j symmetry: 
%iE = 0, 



(17) 



fd'^x I (o"— - co'!— - V"' ^ + C/"^ 
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According to the algebraic renormalization procedure [11], 
Ward identities induce constraints on the most general allowed 
counterterm E' at the quantum level. Once E' is found, one 
can check if it can be reabsorbed in the original starting action 
by a suitable renormalization of fields, sources and parame- 
ters, thereby establishing the renormalizability. One can show 
|10] that E' does not depend on the Lagrange multiplier b", 
and that the antighost c" and the /-valued fields (pf , co-', (pf, W" 
can enter only through the combinations 

K'; = k; + a„c" - gr'^u'^is?" ~ gf'X'To" , 

£/"■ = u;;- + d^m"' , y;' = y;' + a^cp"' , 

A^^' = A^;"' + a^co"' , M'll = y;' + a^cp"' . ( 1 8) 

Imposing the constraints, the most general counterterm yields, 

E' = ooSym 



+fli / d'^x A" 



55 



YM 



" 5A^ 



K"d.c" + V"'m"' - U"'N"J (19) 



V^ 



f IJ f 7J 



with flo, fli two arbitrary parameters. It then turns out that YF 
can be reabsorbed into the starting action (|5]l by a multipHca- 
tive renormahzation ylll^l- At the end, we give the sources 
the following physical values 



<lpM.= ^/iv|,,„,=f§''^§.v, 



u, 



u\u =K\i =K\, ='Ll«/m=0, (20) 
t' \phys f \ phys h'iphys ipliyi ' ^ •" 



in order to recover the physical action (|2]i. 



a forthcoming larger paper flZ|, since the aim of this letter 
is merely to illustrate the relevance of the introduced mass 
operator 

Summarizing, the BRST invariant mass operator cpcp — coco 
fits quite naturally into the theory: it is renormalizable to 
all orders and, moreover, it does not introduce any new 
renormalization constants into the theory. 



IV. THE MODIFIED GLUON AND GHOST PROPAGATOR. 



III. INCLUSION OF A NEW DYNAMICAL EFFECT. 



In a sense, the fields ( (p^'',(p^'',co^' ,co^'' j inti-oduced to local- 
ize the horizon function appearing in ([l]), will correspond to 
the nonlocal dynamics. Once these fields are present, they will 
quite evidently develop their own dynamics at the quantum 
level, which might include further nonperturbative effects, not 
yet accounted for These effects can induce important addi- 
tional changes in the infrared region. More precisely, look- 
ing at the Acp-coupling present at tree level in the action (|2|, 
a nontrivial effect in the cp-sector will get immediately trans- 
lated into the gluon sector. We shall now explore the effects 
of a dynamical mass generation for the (p-fields. This can be 
done by introducing the local composite operator 99 into the 
action (|2]l. Since the horizon condition is in fact equivalent 
with giving a particular value to a dimension 2 Acp-condensate 
i], more precisely fe/"'"A«((p^'' + (p|;')) = -2d {N^ - l)f, 
it does seem to be reasonably fair to consider a possible (pcp- 
condensation. Remarkably, it turns out that this is possible 
while preserving the renormalizability and BRST invariance. 
In order to do so, we try to enlarge the action E by adding a 
massive term of the form f<^<^", with J a new source. First of 
all, for renormalization purposes, we have to add this term in a 
BRST invariant way. Secondly, in analogy with jlOll . we will 
also need a term oc 7-, indispensable to kill potential novel di- 
vergences 0= J^ in the generating functional. We thus consider 
the following extended action: 



2- — 2- + Oqxp? 



•J (pep 



d*x 



d*x 



s{~m^ifi) + p- 



-y((pJ'(pf-co>n + Py 



(21) 



(22) 



with p a new dimensionless quantity and J a new source in- 
variant under the BRST transformation, sJ — 0. We underline 
that the final mass operator, cpcp — coco, is BRST invariant. 
Now, it can be nicely checked that all the Ward identities of 
the previous section are maintained. The final output is that 
the new action ( |2TI ) enjoys multiplicative renormalizability 

An interesting feature is that the anomalous dimension of the 
mass J is not an independent quantity, as it is related to the 
running of the gauge coupling and of the gluon field 111211 . 
We mention already that we will be able to prove that this 
new parameter p is in fact redundant. We postpone this to 



Finally, we come to the main purpose of this letter. We shall 
have a look at the effect on the propagators in the presence of 
the mass operator Ipcp — coco. 



A. The gluon propagator. 

In order to calculate the gluon propagator we only need the 
quadratic part of the action E'. We also replace the source J 
with the more conventional mass notation M^, so that 



e;, = 



|d^x(i (ax - d.A"f + ^ (a,A;)^ +<3>; 



ab 



-I'gif 



abc \(-ifpbc 



2-7^ab ,.ab 



f''''Aim')-M'%y' 



where the limit a ^ is understood in order to recover 
the Landau gauge. The ". . ." stand for the constant term 
~d[N^ ^ 1)Y* and other terms in the ghost and CO, CO fields ir- 
relevant for the calculation of the gluon propagator. Next, we 
integrate out cp and cp, yielding 



Ef> = 



Kab 



d^x-A"A"'f,A: 



ab I 



32 -M2 



^pv 



- d^A 



ti'ib 



Taking the inverse of Af|y and converting it into momentum 
space, we find the following gluon propagator 



<(P) 



where 2'^v(p) — 5, 



-M"- 



p'^+M^p^+2g^Nf 



T.yipW , (23) 



2)(p) 



'pV 



Pl-lPv 



From this expression we make 

three observations: (1) 'D{p) enjoys infrared suppression, 
(2) 2) [p) displays a positivity violation, (3) T) (0) °^ M^, so 
the gluon propagator does not vanish at the origin, which is 
clearly a different result due to the novel mass term propor- 
tional to (pep — coco. 



B. The ghost propagator. 

Let us now turn to the ghost propagator. We first derive the 
gap equation for the Gribov parameter y, useful for the calcu- 
lation of this propagator at one loop. The part of the one loop 




The last expression reveals that the ghost propagator will not 
be enhanced at k^ — 0. Indeed, if we expand (1 — a) in the 
region around k^ = 0, we see 



FIG. 1 : The one loop corrected ghost propagator. 

effective action P''' relevant for this gap equation reads 

4^' = -d{N^-l)j'^ 



2 

(1) 



.-../^K'*-'«) 



Minimizing Fy with respect to y and setting d = 4 leads to 
the following gap equation: 



dV 



3g^N 



J (2% 



1 



{2kY p'^+M^p^ + lg^Nf 



(24) 



We are now ready to compute the ghost propagator at one loop 
order, as depicted in FIG. 1 . The corresponding analytical rep- 
resentation reads 



with 



A^^' 



g"\k) 



kfjkx f d'^q 1 



safe 



1 1 



q^+M^ 



k^ J {2n)'^ (k-q^q'^+M^q^ + lg^Nf 



(25) 



^Ml) 



If we take a closer look at the integral appearing in o, we can 
invoke the gap equation (l24l) in order to simplify (1 — a). By 
virtue of 



Ng- 



2 kfi^v f d'^q 



1 



k^ J (27t)4/+MV + 2g^A^Y* 



J^v- 



q/jqv 



q 



^2^3 A_A 1 



we find 



l-o = Ng 



2kfjk\, 



d^q 

k^ J (27t)4 



ra 



k^-2k-q-M 



21 



' q'^ + M^q^+2g^-N'f 



{k~qf 
^pviq) 



(26) 



1-a 



-Ng- 



2 k^,ky 

k^ 



d'^q -M2 



^Aq) 



(27t)4 q^ q'^+M^-q^ + 2g^Nf 



0{k') 



In 



= — g^M^ 
4« 






8%^ 2y/M'^~8g^N'f 



oik') 



(27) 



We conclude that the ghost propagator keeps displaying a 
1 /k' behavior for k' w 0. It becomes apparent now that, with- 
out the introduction of the new BRST invariant mass term, the 
Gribov-Zwanziger approach would predict a 1/k'^ instead of a 
1 /k' behavior 

V. CONCLUSION. 



In this letter, we have pointed out that the new lattice data for 
the gluon and ghost propagator have a simple understanding 
within the Gribov-Zwanziger approach. This is due to the in- 
troduction of a multiplicatively renormalizable BRST invari- 
ant mass operator 99 — coco, which fits into the theory in a 
very natural way. We hope that the theoretical framework pre- 
sented here will stimulate further investigations, allowing a 
deeper understanding of the propagators in the low momen- 
tum region. We end by noticing that it would be interesting 
to find out what the analogous effects in the maximal Abelian 
gauge might be and also compare those with available lattice 
data liBJ. 
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